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Rotation and reflection symmetries impose that out-of-plane (flexural) phonons of free-standing 
graphene membranes have a quadratic dispersion at long wavelength and can be excited by charge 
carriers in pairs only. As a result, we find that fiexural phonons dominate the phonon contribution 
to the resistivity p below a crossover temperature where we obtain an anomalous temperature 
dependence p oc T^f^ In T. The logarithmic factor arises from renormalizations of the fiexural phonon 
dispersion due to coupling between bending and stretching degrees of freedom of the membrane. 
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Introduction. — The experimental realization of mono- 
layers of graphite, termed graphene, has opened new hori- 
zons in the physics of two-dimensional electron systems 
(2DES) Unlike conventional 2DESj^the low-energy 

electronic bandstructurc of graphene [^, Q, is described 
by a masslcss Dirac equation with velocity v. The 4x4 
matrix structure of the Dirac equation reflects the two 
sublattices of the graphene honeycomb lattice in combi- 
nation with a valley degeneracy due to the presence of 
two Dirac cones within the Brillouin zone. Pioneering 
experiments on this novel two-dimensional electron sys- 
tem have shown that the Dirac nature of carriers induces 
an anomalous integer quantum Hall effect as well as a 
finite conductivity at vanishing carrier density [l|,[2|. 

Recently, it has become possible to experiment on free- 
standing graphene sheets [1, 0| which provide a realiza- 
tion of a two-dimensional (2d) solid. Studies of the sta- 
bility of 2d solids against thermal fluctuations date back 
to early work by Peierls and Landau [1, [§] who pointed 
out the absence of true long-range translational order. 
Much later, it was understood [lO| that nevertheless, 
quasi-long-range translational order can persist up to a 
finite-temperature Kosterlitz-Thouless transition. While 
these works focus on in-plane distortions of the lattice, 
free-standing membranes also support out-of-plane dis- 
tortions. It is believed that there is a low-temperature 
flat phase even in the presence of out-of-plane distortions 
ll|, which gives way to a crumpled phase at high tem- 
peratures [l2| . 

Within the low-temperature flat phase, long- 
wavelength elastic distortions, both in-plane and 
out-of-plane, can be described by the appropriate elastic 
Lagrangian density 11 , 3] , 



(1) 



in terms of the mass density pq , the out-of-plane distor- 
tions h{r), and the strain tensor Uij = ^[diUj + djUi + 
{dih){djh)]. Here, u(r) denotes the in-plane distortions. 
The elastic constants A and p characterize the in-plane 
rigidity of the lattice, k the bending rigidity. Both the 
absence of a (V/i)^ term in the elastic Lagrangian and 



the appearance of the out-of-plane distortions h{r) in the 
strain tensor are direct consequences of the rotational 
symmetry of the membrane in the embedding space. 

The elastic Lagrangian in Eq. H]) encapsulates a dis- 
tinct difference between in-plane and out-of-plane (flex- 
ural) phonons. Indeed, to quadratic order in the dis- 
placement fields h{r) and u(r), both longitudinal and 
transverse in-plane phonons have a linear dispersion 



,(0 



= w^'-'q and Wq = w^'^g with group velocities 

— [{2p + X) / po]^^^ and w^*^ = [p/po]^^^. In contrast, 

flexural phonons obey a quadratic dispersion tjq'*'' = aq^ 
1/2 

with a = [k/po] , which is a consequence of rotational 
symmetry. In-plane and flexural phonons also differ in 
their coupling to the charge carriers of graphene mem- 
branes. While the coupling is conventional for in-plane 
phonons, the reflection symmetry h — > —h demands that 
out-of-plane displacements enter only quadratically into 
the Dirac Hamiltonian. Consequently, charge carriers 
can excite flexural phonons only in pairs. 

Due to these differences, we find that flexural phonons 
dominate the phonon contribution to the resistivity of 
free-standing graphene membranes below a crossover 
temperature T^. Indeed, the transport scattering rate 
of Dirac fermions diverges logarithmically for a strictly 
quadratic dispersion of flexural phonons. This divergence 
is cut off by the coupling between bending and stretching 
degrees of freedom of the membrane, as captured by the 
elastic Lagrangian Eq. ([T]). At finite temperature, this 
coupling renormalizes the bending rigidity of the mem- 
brane, inducing a stiffening of the flexural-mode disper- 
sion at long wavelengths. Including this physics within a 
simple one-loop RG, we find that the contribution of fiex- 
ural phonons to the resistivity of graphene membranes 
scales as T^/^ InT. 

Graphene. — The bandstructurc of graphene is well ap- 
proximated by the tight-binding Hamiltonian 



H=-tj: 
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on a honeycomb lattice. Here, t is the hopping matrix 
element, Ci annihilates an electron on lattice site i, and 
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only nearest-neighbor hopping has been included. The 
2d hexagonal lattice consists of two identical sublattices 
A and B, and thus two sites per unit cell. We denote 
the vectors connecting a B site with the neighboring A 
sites as = a (—1,0), 63 = a (1/2, •\/3/2) and 63 = 
a(l/2,— v^/2) where a is the bond length. The band 
structure of the Hamiltonian in Eq. ([2]) has zero energy 
(corresponding to the Fermi energy at half filling) at two 
inequivalent points in the Brillouin zone, which we choose 
to be at k_(_ = ik^,, with k^, = 27r/(3-\/3a) (-y/S, 1). In 
the vicinity of these Dirac points, the spectrum is de- 
scribed by the 4x4 Dirac Hamiltonian 



H = hv'E-k 



(3) 



with velocity v = 3ta/2. The 2d wavenumber k is mea- 
sured from the Dirac point. The Hamiltonian in Eq. ^ 
acts on four-component spinors {u\ ^, j., ^, j.) 
of Bloch amplitudes in the space spanned by the sub- 
lattices (A/B) and Dirac points (+/—). The matrices 
'^x.y = (g) a,j. y denote components of a vector S. 
(Hj and tr^ are Pauli matrices acting in the spaces of 
the Dirac points and the sublattices, respectively). It 
is also useful to introduce a corresponding vector A by 
A^ 



-x,y - ^x,y cTz and A^ = H^ 
Electron-phonon coupling.— 



) and = il^ (X) ctq 

-The dominant electron- 
phonon coupling arises from distortion-induced modifica- 
tions of the bond lengths and hence the hopping ampli- 
tude [15, H, 17|. The electron-phonon coupling can then 
be described in terms of a fictitious gauge field A(r) en- 
tering into the Dirac Hamiltonian, H = z;S-(p-f eA(r, t)), 
where p denotes the momentum. Expressed in terms of 
the strain tensor Uij, the gauge field A(r) takes the form 

M 



eA(r,i) = H, 



t da 



2 {^xx '^y 
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The combination of components of the strain tensor en- 
tering into A(r) reflects the symmetry of the underlying 
honeycomb lattice. The factor of H^ implies that the as- 
sociated flctitious magnetic field points in opposite direc- 
tions at the two Dirac points, as required by time-reversal 
symmetry. From the definition of the strain tensor, we 
directly conclude that electrons couple linearly to the in- 
plane phonons and quadratically to flexural ones. The 
corresponding electron-phonon vertices are depicted in 
Fig. [11 

At sufficiently low temperatures, only long-wavelength 
phonons contribute to the resistivity and scattering be- 
tween different Dirac cones can be neglected. Thus, we 
can restrict attention to the vicinity of, say, the Dirac 
point k_|_ where the Hamiltonian reduces to a 2 x 2 Hamil- 
tonian in the A — B sublattice space. 

The electron-phonon coupling H^^ corresponding to 
Eq. can be expressed in second quantization af- 
ter expanding the in-planc and out-of-plane distortions 
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FIG. 1: Electron-phonon vertices, a) Coupling of electrons 
to in-plane phonons. b) Coupling of electrons to flexu- 
ral phonons. Straight lines correspond to electrons, wavy 
(dashed) lines to in-plane (flexural) phonons. 



into Fourier series as u(r) = X^q^qe''^ '' and h{r) = 
J2q e*'* '", decomposing Uq into longitudinal and trans- 
verse components. 



u^*^z X q 



and quantizing the amplitude of the distortions as 



2Muj\ 
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Here M is the atomic mass and a'c^^ the annihilation oper- 
ator of a phonon of type v [v = l,t, h) with wavenumber 
q. This yields 

^ep- E E Ki'^^4,kCA,k'(<+«^'^^)5k',k-q 
k,k',q y—l.t 



+ E ^q,q''^kk'^A,k' ( 

i,q' 



k,k' ,q,q' 

X 
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k-q-q' 



h.c. 



in terms of the coupling terms 
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2Muj. 
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q^q 
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where e = {3a/A)dt/da, ^ denotes the creation oper- 
ator of an electron in a Bloch state in sublattice B and 
wave-vector k, and </> {(/)') is the angle of q (q') with 
respect to the x-axis. 

These considerations give rise to an interesting compe- 
tition: In-plane phonons are strongly coupled to electrons 
(the vertex is first order in the phonons) but they have a 
linear dispersion and hence a linearly vanishing density of 
states at small energy. In contrast, flexural phonons are 
weakly coupled to electrons (their vertex is second order 
in the phonons) but their dispersion is quadratic with 
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constant density of states. In the following, we quanti- 
tatively analyze the consequences of this competition for 
the temperature dependence of the resistivity. The latter 
is determined by 



(9) 



{lyp is the electronic density of states at the Fermi level) 
in terms of the transport scattering rate 

— = \Mff {1- cose) S{Ef ^ E,±huj) . (10) 



of Dirac fermions due to absorption (or emission) of 
phonons. Here, Q is the scattering angle of the Dirac 
fermions, \i) and |/) the initial and final scattering 
states, Ei and Ef the initial and final electronic ener- 
gies, hio the energy of the absorbed (emitted) phonons. 



and M 



op I 



We consider the doped regime with Ep ^ fcgT and 
employ the conventional quasielastic approximation of 
neglecting the phononic contribution to energy conserva- 
tion, as the typical phonon energies are small compared 
to Ep. Moreover, we focus on sufficiently low tempera- 
tures so that we can restrict attention to the quadratic 
(linear) region of the dispersion of flexural (in-plane) 
phonons. 

Flexural phonons. — Scattering of Dirac fermions by 
flexural phonons requires absorption (or emission) of 
two phonons, say with wavcnumbcrs q and q'. Thus, 
the corresponding initial and final scattering states are 
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k,(T) (g) \ni^^) and |/) = |k + q + q\a' 
— 1). where nf*-* 



1, n 



denotes the Bose distribution function of the flexural 



phonons and |k, cr) = l/y/2(ae '^cj. 



vac) IS 



a Dirac fermion state state with momentum k and chi- 
rality cr. (Here, ^ is the angle between k and the x- 
axis and |vac) is the electronic vacuum). Thus, = 

iK^.q' P [1 " f^f^' cos(C + e' + 20 + 20')] /2. The summa- 
tion over final states requires integration over q and q', 
while averaging over the direction of the incoming elec- 
tron will suppress the oscillatory term in |Af^jp. Using 

ai^^ln^^^) = {ni^^)'/'\ni^^ -I), E, 
Ej = hv\\s. + q + q'l, we obtain 



hvk — Ep, and 



Jh) 



(h) 12 



(l-cos0)nWn(^)<5($q+qO, (11) 



Ep and \V^^,\^ 



where $q+q' — ?iu|k -|- q -I- q' 

?i^e^/(16Af^a^) is independent of wavenumbers due to 
the quadratic dispersion of flexural phonons. The rate 
in Eq. pip is formally singular at small q, q' since both 
Bose distributions diverge as T jq^. Rcscaling momenta 



FIG. 2: One-loop correction to the bending rigidity due to 
the effective interaction between flexural modes. 



by vT and introducing a cutoff q^ at small wavenumbers 
(to be specified below), we obtain the scattering rate 



1 C^''' ^'^k 



(h) 



5/2 
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huj 



In 



k^T 



huj, 



(h) 

9c , 
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where C^"^) = e^hkj./{4poEpnu}[''^). The unusual I 
scaling implies that scattering from flexural phonons 
dominates the phonon contribution to the resistivity at 
low temperatures. Indeed, the conventional phonon con- 
tribution to the resistivity due to in-plane phonons scales 
as T^, which is the direct two-dimensional analog of the 
law in bulk solids. 

Eq. (fT2|) shows that the scattering rate from flexural 
phonons diverges logarithmically for a strictly quadratic 
phonon dispersion. For clean elastic membranes, a low- 
momentum cutoff arises from the coupling terms between 
bending and stretching degrees of freedom of the mem- 
brane which are included in the Lagrangian Eq. p]). 
These lead to long-wavelength corrections to the elastic 
constants and hence the phonon dispersions ll| . Indeed, 
it is these renormalizations which are responsible for the 
stability of the flat phase at low temperatures. 

In order to study the renormalization of the bend- 
ing rigidity by the coupling of bending and stretching 
modes, we integrate out the in-plane distortions in Eq. 
(H]) and obtain an effective energy functional for the flex- 



ural modes alone [ll| 



F = 



dr 



(13) 

where = 4/i(^ + + A) and P^^ = 

(^q/ ^/3fe ^/ ^fc)/^^ ^ transverse projector. This effec- 
tive energy functional contains a four-leg interaction be- 
tween flexural distortions. Treating this quartic term 
to one-loop order (depicted in Fig. [5]) in a conventional 
momentum- shell RG, one obtains the flow equation 



dK 
dq 



3 
16^ 



nq^ 



(14) 



with q denoting the running shell wavevector. 

The flow equation Eq. is readily solved, and the 
resulting scaling of k with q is 



K{q) = K^/lTq^J(f 



(15) 
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in terms of the temperature-dependent momentum scale 
= [iKQk^T/{d,TTK^)Y/'^. Thus, the flexural phonon 
dispersion is quadratic for q ^ q^. In contrast, for q <C 
q^., thermal fluctuations effectively stiffen the membrane 
and, within the simple one-loop analysis presented here, 
we find a renormalized dispersion uj^j^q = [Kq^)^/'^q^/'^ . 
It is this renormalization of the flexural-phonon disper- 
sion at long wavelength which removes the singularity in 
the phonon scattering rate Eq. (|12p. We therefore iden- 
tify the low-momentum cutoff qc entering into Eq. (jl2p 
with this momentum scale. 

As a result, we find that the temperature cancels from 
the argument of the logarithm in Eq. (fT2|) . This anal- 
ysis is approximate in that it is restricted to one-loop 
order and that it neglects the flow of the elastic con- 
stants for stretching deformations. While the stiffening 
against bending deformation would survive inclusion in 
a more refined treatment, the precise power-law depen- 
dence of 5c on temperature would change. As a result, 
the argument of the logarithm would become tempera- 
ture dependent, resulting in an overall T^/^lnT scaling 
of the flexural phonon contribution to the resistivity. 

Crossover temperature. — In order to estimate the 
crossover temperature Tx below which phonon scatter- 
ing of Dirac fermions is dominated by flexural modes, we 
note that an analogous calculation yields 



huj 



(16) 



for the transport scattering rate from in-plane phonons 
[v = l,t). Thus, the crossover temperature obtained from 
a comparison of Eqs. ^2]) and becomes 



1 



In ( k^T /huj, 
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huj 



(t)5 



2/3 



K') 
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(17) 



This result is independent of the Fermi energy (i.e. the 
doping level) and with typical parameters for graphene 
[2lj . we obtain ~ 70 K. There is therefore a significant 
temperature range accessible in experiment over which 
our predictions for the resistivity can be tested. 

Conclusions. — For clean graphene membranes, long- 
wavelength renormalizations of their elastic properties 
due to thermal fluctuations are crucial in order to ob- 
tain a finite transport scattering rate and hence resistiv- 
ity p (X T^/^lnT from scattering by flexural phonons. 
In the presence of disorder, the elastic properties of the 
membrane are renormalized even at zero temperature 
Til . 20| . Indeed, the recent experimental observations of 
rippling Q suggest that disorder, inside or close to the 2d 



membrane, exists in present graphene membranes. Such 
disorder-induced renormalizations of the elastic moduli 
may compete with the renormalizations by thermal fluc- 
tuations and lead to a temperature-independent satura- 
tion of the cutoff qc at low temperatures. A detailed 
study of the effects of disorder in this context remains an 
important topic for future research. 
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